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ON AN APPROACH TO CONSTRUCTING STATIC BALL MODELS 

IN GENERAL RELATIVITY 1 
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Abstract - An approach to construction of static models is demonstrated for a fluid ball. Five examples are 
considered. Two of them are exact solutions of the Einstein equations; the other three are connected with the Airy 
special functions, the hypergeometric functions and the Heun functions. 
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One of the most important problems in general relativ- 
ity is that of finding exact solutions of the gravitational 
equations. Unfortunately, this task is not easy. An- 
other way is an analytic construction of solutions to the 
gravitational equations with a certain physical interpre- 
tation. 

Further on we take the metric in Bondi's form 



ds 2 = G{r) 2 dt 2 + 2L(r)dtdr 
-r 2 (d0 2 + sin 2 (0)dtp 2 ) 



(1) 



where G 2 ,L are metric functions, r is a radial variable 
and 9, ip are spherical angles; the speed of light and 
Newton's constant of gravity are put equal to unit. 

The gravitational field is described by the metric ten- 
sor gn~, which can be found from Einstein's equations 



Rik — -^Rgik — —xTik, 



(2) 



where k — 0,1,2,3; Rik is the Ricci tensor, R is 
the scalar curvature of space-time; x is Einstein's grav- 
itational constant. 

The energy-momentum tensor (EMT) of Pascal's 
perfect fluid can be written as 



T ik = (p (r) + ju (r)) • UiUk - P (r) g ik , 



(3) 



where p{r) is the pressure, fjb[r) is the mass density, it, 
is the 4-velocity. 
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The gravitational equations in dimensionless vari- 
ables can be reduced after elementary transformations 
to the form 



e(x) = 1 • / n(x) ■ x 2 dx; 
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(4) 
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L-4 ) r; = o: 



p 



Xxp+ ) • (fi + p), (6) 
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where x — r/R is the dimensionless radius; < x < 1; 
differentiation in x is denoted by a prime; R is the 
radius of the astrophysical object; x = xr • i? 2 ; 



e{x) 



G 2 (x) 
L 2 (x)' 



(7) 



2. Reduction of the Einstein 
equations 

Eq (5) can be reduced to an oscillatory-type equation 
d 2 G 



where 



fK ., , n 2 (((y))G = 0, 
n 2 = -d(<P/y)/dy, 



y = x 2 , C is a new variable: 



d( = ydy/{2yfe). 



(8) 

(9) 
(10) 
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Now we introduce e = 1 — $ as in [T] , where 

$ = (x/(2Vv) J vVydy, (n) 

where \ — xR 2 , $ is an analog of Newton's gravita- 
tional potential. 

Further we present the function O 2 in the form of a 
power series: 



(12) 



Now we can find both <I> and fj, in the general case 
from (8) and (10) with help of (11) 



$ = ( Mo /3) y -^^- y » 
i n + 1 



(13) 



n=0 



M / Mo = i-(i/Mo)E^rr a ^ n+1 ' ( 14 ) 

* — ' 71+1 

71=0 

where Mo is the central mass density and the coefficients 
a n are constants to be found from boundary conditions. 

3. Construction of solutions 

to the gravitational equations 

Consider the power series (11). At first we will take all 
coefficients a n equal to zero. After that we will consider 
only ao 7^ 0, further ag =^ 0, a\ ^ and so on. In each 
case we will construct the corresponding mass density 
and make an attempt to find the function G from (7). 
Here we must remark that the metric functions </oo = G 2 
and g w = g 01 = L = G/^/e. 

1. If we take all a n = 0, then O 2 = and Eq.(8) is 
transformed to 



fG 



= 0. 



(15) 



The solution is G(C) = Ci • C + C 2 , where C\ , C2 are 
integration constants. From (14) we have u = Mo- In 
other words, it is the Schwarzschild interior ball model. 

2. If we include only ao ^ 0, then O 2 = 0§ = const 
and 

G oc cos(flo( + ipo). (16) 



In this case 



M = /jq - 3a y = /Jo - 3a x 



(17) 



i.e. we have a parabolic distribution of the mass density 
(see H). 

We must say that in these cases the approximate 
solutions coincide with the exact well-known solutions 
of the Einstein equations. 

3. Now we take £ approximately as £ ps y/2 (C(0) = 
0), because there is the difficulty in determining the 



variable C via y. Here we have an approximate estimate 
because the variable is y « 1. In this case the origins 
of the two variables y and £ are glued. 
The further approximation will be 



fi = a + aiy ps a + 2ai£, 



and Eq.(7) can be written as 
d 2 G 



rfC 



2 +( ao + 2aiC)G = 0, 



(18) 



(19) 



while the mass density is 

(j, ps /to - 6a C - lOaiC 2 

« Mo - 3a a; 2 - (5/2)aix 4 . (20) 
The solution of Eq.(19) is 



G = CiAiryAi 
+C 2 AiryBi 



.fao + 2ai£ 



(2ai) 2 / 3 

a + 2aiC 
(2ai) 2 / 3 y 



(21) 



where AiryAi and AiryBi are the Airy special func- 
tions, and Ci,C2 are constants. 

4. The following real solution with 



and 



fr pa a + 2aiC - 4a 2 C 



/j ps /jo - 6a C - lOaiC 2 + (56/3)a 2 C 3 



(22) 



Mo - 3a x 2 - (5/2)a lX 4 + {7/3)a 2 x 6 (23) 



can be written as the linear combination of hypergeo- 
metric functions: 

G = [C 1 hypergeom(ax, /3i,7(C))+ 

C 2 hypergeom(a 2 , /3 2 , l(C))(^ a 2C - 

xe^^(C), (24) 
where C\,C 2 are constants, 

«i = — (ai/4 + a a 2 - 2a 3 / 2 )/ (8a 3 / 2 ); 

/?! = 1/2; a 2 = ai + 1/2; 

ft = ft + 1; 7(0 = (4a 2 C - ai) 2 /(8^ /2 ) ; 

«(C) = C(ai/2-aaC)/VS. 
5. One more real solution can be found for 

fi 2 wa + 2aiC-4a 2 C 2 -16a 4 C 4 (25) 

M ps /jo - 6a C ~ lOaiC 2 + (56/3)a 2 C 3 

-(352/5)a 4 C 5 « Mo - 3a a: 2 - (5/2)aia; 4 

+ (7/3)a 2 x 6 + (ll/5)a 4 ir 10 . 

This solution is written as a linear combination of 
HcunT functions 
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G = dHeunT(a, -/3, 7, -b-Q- exp(ip(())+ 

+C 2 HeunT{a, /3, 7, b ■ £) ■ exp(-7p(Q), (26) 
where C\,Ci are constants, 

a = 3 2 / 3 (4a 4 ao + 4) / '(ma 4 / 3 ); /3 = 3^/(4^1), 

7 = 3 1 / 3 a 2 /(2af ),6 = (2/3)3 2 / 3 a y 6 , 
V-(C) =C-(3a 2 + 8a 4 C 2 )/(6V^). 

The pressure can be found from the equation Eq.(6) 
and the metric function L from (7). 

4. Summary 

In conclusion, we must note that an approach to the 
construction of the ball static models is demonstrated 
in this paper. This approach is based on the reduction 
of gravitational equations to the oscillatory-type equa- 
tion and the using the expansion in the power series 
the function which plays a role of the frequency. Main 
difficulty is to find the new variable through the dimcn- 
sionless radial variable. Five examples are considered. 
Two from them are the exact solutions of the Einstein 
equations for a fluid ball. The third, fourth and fifth 
examples are connected with the special Airy functions, 
with the hypergeometric functions and the HeunT func- 
tions. 
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